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The WDVV equations of associativity in 2-d topological field theory are 
completely integrable third order Monge-Ampere equations which admit bi- 
Hamiltonian structure. The time variable plays a distinguished role in the 
discussion of Hamiltonian structure whereas in the theory of WDVV equa- 
tions none of the independent variables merits such a distinction. WDVV 
equations admit very different alternative Hamiltonian structures under dif- 
ferent possible choices of the time variable but all these various Hamiltonian 
formulations can be brought together in the framework of the covariant the- 
ory of symplectic structure. They can be identified as different components of 
the covariant Witten-Zuckerman symplectic 2-form current density where a 
variational formulation of the WDVV equation that leads to the Hamiltonian 
operator through the Dirac bracket is available. 



1 Equations of associativity 

In 2-d topological field theory Witten || has shown that model indepen- 
dent n-point correlation functions follow recursively from the 2- and 3-point 
correlation functions that serve to define a non-degenerate flat metric and 
structure functions of a Frobenius algebra. These are the principal objects 
in this theory and they can be expressed as third derivatives of a generating 
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function 

d 3 F 

where F is the free energy ||. One of the independent variables, say t 1 , is 
singled out to serve in the definition of the metric 

C Uj = T)ij (2) 

which is assumed to be non-degenerate. The structure functions satisfy a 
Frobenius algebra and the conditions of associativity 

Cmi[j C k ]ln = (3) 

result in third order Monge-Ampere equations which are known as WDVV 
equations. Indices enclosed by square parantheses are skew-symmetrized. 

Dubrovin |4[] has given a systematic account of WDVV equations of asso- 
ciativity which are completely integrable systems. They admit bi-Hamilton- 
ian structure H , which also provides proof of their complete integrability 
through the theorem of Magri 0. However, in the discussion of Hamiltonian 
structure time, which can be identifmed to be any one of t\i ^ 1, is neces- 
sarily singled out whereas in the general theory of WDVV equations no such 
distinction exists. To give an example, consider the free energy 

Fl = l -^fe+ l -t\ef + f{t\e) (4) 

that through the identification 

t 2 = x } t 3 = t (5) 
results in Dubrovin's equation of associativity 

fttt + fxxx fttx — ftxx 2 = 0- (6) 

By a trivial interchange in the roles of t and x, alternatively if in place of (|^) 
we make the identification 

t 2 =t, t 3 = x, (7) 

then we are led to 

fttt ftxx — fttx 2 + fxxx = (8) 
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which should be the same equation of associativity as there is no distinction 
between the independent variables t and x that stemms from 2-d topological 
field theory itself. However, from the point of view of Hamiltonian struc- 
ture Eqs.(|6|) and @ are radically different. The bi-Hamiltonian structure of 
Eq. (@) which is based on the results of || and || was given in [||] and in this 
paper we shall present the bi-Hamiltonian structure of Eq.(|8|). We shall find 
that the Hamiltonian operators appropriate to eq.(||) are quite different from 
those obtained for eq.(|6|), however, we shall further show that they can be 
identified as different aspects of the same structure when we consider them in 
the framework of the covariant Witten-Zuckerman [|T0| formulation of sym- 



plectic structure. Considerations of covariant symplectic structure take the 
variational formulation as their starting point. For WDVV equations only the 
Lagrangian that yields their second Hamiltonian structure through Dirac's 
theory of constraints []TTj| is known and therefore this part of our discussion 



will necessarily be restricted. We shall show that the second Hamiltonian 
operators for eq. (P) and (|]) are simply the inverse of the t and x components 
of the Witten-Zuckerman symplectic current 2-form for only one of these 
equations. 

Exactly the same situation holds for the WDVV equation that follows 
from the free energy 

f 2 = \{t"f e + V (t 2 ) 2 - Ipf e + f{t\ t 3 ) (9) 

which through the identification leads to 

fttt fttt fxxx fttx ftxx fttt ftxx fxtt fxxx fxtt fxxt (10) 

and the bi-Hamiltonian structure of this equation was presented in ||. With 
the same free energy the opposite identification ([!]) results in 

fttt fxxx fttx ftxx fxxx fxtt fxxt fttt ftxx fxtt fxxx (11) 

as the WDVV equation of associativity. Once again we shall find that the bi- 
Hamiltonian structure of eq.flTTD is different from the earlier results obtained 
for eq. QTOl) but they can be recognized as different components of the closed 
conserved current 2-form in the covariant formulation of symplectic structure. 
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2 System of evolution equations 



In earlier literature |L^] Monge-Ampere equations which are second order 
partial differential equations consisting of an appropriate sum of linear terms 
and the Hessian were discussed in the framework of linear equations proper. 
The reason for this lies in the fact that the initial value problem for hyper- 
bolic Monge-Ampere equations is identical to that of the second order linear 
equation [13] in spite of the severe non-linearities introduced by the Hessian. 



WDVV equations of associativity are of third order but they consist of a sum 
of linear terms and the minor determinants of the third order Hankelian 

H(u,t,x)=detminor\ fm fttx ftxx ) (12) 

I Jttx Jtxx 1 XXX J 

closely analogous to the case of Monge-Ampere equations. Therefore it seems 
reasonable to conjecture that the initial value problem for WDVV equations 
is qualitatively the same as the third order linear equation. This is partic- 
ularly important in the discussion of the Hamiltonian structure of WDVV 
equations where we shall use the techniques originally developed for real 
Monge-Ampere equations 0]. For this purpose we need to cast the WDVV 
equation into the form of a triplet of evolution equations. Introducing the 
usual auxiliary variables 

a = fxxxi b = fxxti C = fxtti (13) 

we have in place of the WDVV equation (||) the set of evolution equations 



(14) 



a t = 


b x , 


b t = 


Cx ; 


Ct = 










c 2 — a 


e = 


b 



which consist of equations of hydrodynamic type |y| • We note that such a 
decomposition of eq.(|8|) is not unique but this particular choice of auxiliary 
variables is useful because the Hamiltonian structure of equations of hydro- 
dynamic type is a well developed subject. In the case of eqs.(|i"4"|) the system 
is linearly degenerate and it will be necessary to use the results of Ferapontov 
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on the Hamiltonian structure of non-diagonizable equations of hydrody- 
namic type for which Riemann invariants do not exist. From eqs.(|l4]) it is 
manifest that a, 6, c are conserved densities and there are two others 



V 
£ 



bD~ l c, 
cD^bD-h 



D~ x aD~ x b 



(15) 



which consist of momentum and energy. There exist no further conserved 
quantities of hydrodynamic type as eqs.(|l4l) are non-diagonizable. Here and 
in the following D~ x stands for the inverse of the total derivative operator 



D = d/dx and its precise definition can be found in fl7| . 



3 Bi-Hamiltonian structure 



We shall first state the principal result we shall present about eq.(||). 
Theorem 1: 

Eqs. fllil) can be written bi-Hamiltonian system 

J 5 H = Ji 8Hi 



(16) 



where 8 denotes the variational derivative with respect to a, b, c and the 
Hamiltonian operators given by 



Jn 



/ 2aD + 2Da 3bD + b x 2cD 
3bD + 2b x cD + Dc eD 
2Dc De -3D 



Ji 



V 

/ —D 3 













V 



o 



d\d> d\d^d + dId\dj 



(17) 



(18) 



are compatible so that according to Magri's theorem |7| we have a completely 
integrable system. 

The Hamiltonian functions that yield the equations of motion are integrals 
of the densities TLq = b and Hi = £ respectively. 
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4 Spectral problem 



The general framework for casting WDVV equations into the form of a spec- 
tral problem was given by Dubrovin [|J. In the case of eq.@, or rather 
eqs. (|14D these considerations yield the Lax pair 







zAm, 


% = 


zBV, 




(° 


1 








(00 




A = \ a 





!)■ 


B = 




b 


:) 


\b 











{c 1 





where z is the spectral parameter and the compatibility conditions 

A t = B x , [A,B] = 

are satisfied by virtue of eqs. ( [Ill ) . From the compatibility conditions it follows 
that the roots w 1 ,^ 2 ,^ 3 of the characteristic equation 

det(XI-A) = \ 3 -c\ 2 -aX + ac-b 2 = (20) 

are conserved Hamiltonian densities for eqs.(|H|). 



5 First Hamiltonian structure 

The easiest way for finding the first Hamiltonian structure of eqs.flli]) is to 
transform to a set of new dependent variables which consist of the roots of 
the cubic ( |20"D because the Hamiltonian operator then assumes the form of a 
first order homogeneous operator with constant coefficients. The relationship 
between the roots of this cubic to the hydrodynamic variables a, b, c is given 
by 

a = —j3, a = u 1 +u 2 + u 3 , 

b = =fV7 - a 13, 13 = u l u 2 + u 2 u 3 + u 3 u 1 , (21) 

c = a, 7 = u^u 2 u 3 

according to the formulas of Viete. In these variables the equations of motion 
fllip assume the form 

« = 1,2,3 (22) 



^7 - a/3 
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of Hamilton's equations with the Hamiltonian operator 

j/i -i _i\ 

J =- -1 1 -I D (23) 




and the Hamiltonian density 

Ho = b = ^u l u 2 vfi — (u 1 + u 2 + m 3 )(-u 1 'u 2 + u 2 u 3 + vfiu 1 ). 

A lengthy but straight-forward transformation of this operator into the orig- 
inal auxiliary variables yields the result in eq . (|17D . The verification of the 
Jacobi identities for Hamiltonian operator (|TT| ) is straight-forward. For this 
operator the conserved quantity \a results in the trivial flow while c is a 
Casimir. 



6 Variational principle 



The most direct way of obtaining the second Hamiltonian structure of eq. (H) 
is through the construction of a variational principle. The auxiliary variables 
a, b, c are not best suited for this purpose, instead it turns out that in terms 
of 

P = fx, Q = ft, r = f tt , (24) 

the Lagrangian is given by a simple local expression. The equations of motion 
are 

Pt = <lx, 



Qt = r, 

r — r x ~ Pxx 



(25) 



Qxx 

and it can be readily verified that the variational principle with the La- 
grangian density 

£ = ^Px Pt + Qx r x q t - ^q 2 x r t + p q xx + ^r 2 q xx (26) 

yields eqs.(^5]). We note that this Lagrangian is linear in the velocities so that 
the Hessian vanishes identically and we have a degenerate Lagrangian system. 
The passage to its Hamiltonian formulation requires the use of Dirac's theory 
of constraints [ pT|j . 
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7 Dirac bracket 



The Dirac bracket which replaces the Poisson bracket for systems subject 
to constraints plays a central role in the construction of the Hamiltonian 



operator for integrable systems |lq| . This construction is directly applicable 
to the WDVV equation (§) as we shall now detail. Since the Lagrangian ( p6|) 
is degenerate the canonical momenta cannot be inverted for the velocities 
and following Dirac we introduce the definition of the momenta 

01 = 7T p - \p x , 

02 = TTq ~ q X T x , (27) 

03 = ?Tr + 

as primary constraints. Calculating the Poisson bracket of the constraints 
{0i(x),0i(?/)} = \8 x {y - x) - \8 y {x - y) , 

{M x )><f>s(v)} = o, , 2g s 

{02 (x), 02 (y)} = r x 5 x (y - x) -r y 5 y (x - y) , " 
{02 (x), 03 (y)} = q x 5 x (y - x) + q y 5 y (x - y) , 

{h(x),Mv)} = o, 



we find that the constraints (|27j) are second class as they do not vanish mod- 
ulo the constraints themselves. This is the case with almost all completely 
integrable systems. The total Hamiltonian density of Dirac is given by 

H T = H 1 + T* =1 d ^ H x = ~-r 2 q xx -pq xx (29) 

where d are Lagrange multipliers. The expression for TCi is obtained from 
the Lagrangian fl2"E| ) by Legendre transformation. The conditions that the 
constraints are maintained in time {(pi(x), H?} = give rise to no further 
constraints but rather determine the Lagrange multipliers 

2 

1 2 3 T x ~ Vxx 

c =q x , c =r, c = 

q X x 

and we have no secondary constraints. From eq.(p9|) we find that the total 
Hamiltonian is given by 

2 — p xx / 1 \ 1 

Ht = q x n p + r 7i q + — — I vr r + -q 2 J + -p x q x (30) 

q xx \ Z J Z 
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in terms of the full set of canonical variables. For systems subject to sec- 
ond class constraints we can solve the constraints to eliminate the canonical 
momenta because in Dirac's theory second class constraints hold as strong 
equations. As a result the total Hamiltonian density is simply TC% which up 
to a total derivative is the same as S in eq.(|l5l). 

Given any two smooth functionals A, B the Dirac bracket is defined by 

{A(x),B(y)} D = {A(x),B(y)} 

-J{A(x),^)}J lk (z,w){Mw),B(y)}dzdw 

where J lk is the inverse of the matrix of Poisson brackets of the constraints. 
From the definition of the inverse 

{(f)i(x), (f) k (z)}J kj (z,y)dz = 5j5(x-y) 

we end up with a set of differential equations for J tk which can be solved to 
yield 

J 11 (x,y)= -9(x-y), 

J 12 (x,y) = 0, 

J 13 (x,y)= 0, 

J 22 (x,y) = 0, ( 32 ) 

J 23 (x,y) = —^S(x-y), 

J 33 (x,y)= -2^8 x (x -y) + Tx % xx S(x - y) 

Qxx Qxx 

where the Heaviside unit step function is denoted by 9. 

8 Second Hamiltonian structure 

The transition from the Dirac bracket to the Hamiltonian operator is given 
by 

{u\x) lU k (y)} D = -T k (x,y) = -J ik {x)8{x-y) (33) 

since the Poisson brackets of u % vanish and the constraints are linear in the 
momenta. From eqs.(p3|) and (|32|) it follows that the Hamiltonian operator 
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corresponding to the degenerate Lagrangian ( |26| ) is simply 



■h 



(D-* 




V 









J_ 





J_ 

(fes 9cca; 



(34) 



and the proof of the Jacobi identities for the Hamiltonian operator fl34|) fol- 
lows from the fact that according to eq. ( |33"D it is simply a reformulation of the 
Dirac bracket for which there is a general proof of the Jacobi identities [H] . 
The Hamiltonian operator (B3) looks non-local but this is only superficial 



and related to the choice of variables p, q and r. If we revert to the original 
auxiliary variables fllTf ) 



a = p a 



the Hamiltonian operator (34) is transformed to the form (18) which is a local 
homogeneous third order operator of hydrodynamic type that was studied in 



9 Symplectic representation 



The symplectic formulation of eqs. (|25|) that provides the dual description to 



the Hamiltonian operator fl34]) requires the inverse of this operator. However, 
we have just seen that this Hamiltonian operator is derived from the Dirac 
bracket which in turn was obtained from the inverse of the Poisson bracket 
of the constraints. The matrix of symplectic 2-form density can therefore be 
obtained directly from eqs.(P%|) and we get 



DO 

Uij = - | Dr x + r x D -q xx | (35) 
q xx 



which can be verified to be the inverse of (0). The symplectic 2-form density 
is then given by 

lu = ~^dp A dp x - r x dq A dq x + q xx dq A dr (36) 
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and since this 2-form is closed, using Poincare's lemma we can write 

uj = da, , 

(37) 

a= \p x dp-rq xx dq 

in a local neighborhood. In terms of / that enters into the original formula- 
tion of the WDVV equation (|8|) the symplectic 2-form PB| ) reduces to 

u = 2 f tx df tt A df tx - \- df x A df xx (38) 

using the definitions (|24] ) and discarding a total derivative. 



10 Darboux's theorem 



The transformation of the third-order Hamiltonian operator 
canonical form 

/ 1 
Ji = - 1 I D 

I i o o 



H) to the 



(39) 



required by an as yet unproved generalization of Darboux's theorem for third 
order Hamiltonian operators is usually given by a differential substitution 
where the Casimirs of J\ play a central role. For the Hamiltonian operator 
( P~8| ) the Casimirs are 

s 1 = D^b, 



bD 



(40) 



and the differential substitution required by Darboux's theorem is obtained 
by inverting eqs.(||]). Indeed it can be readily verified that under this trans- 



formation of variables the second Hamiltonian operator (18) is transformed 
to the form (|39|) which is proof of Darboux's theorem for this case. Finally we 
note that in terms of the Casimir variables the equations of motion assume 
the form 



(41) 
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of an integrable coupled dispersive system. 



11 Summary of earlier results 



In order to compare the results on bi-Hamiltonian formulations of eqs.@ 
and (P) we need to summarize the results of |J on eq.([|). In this case 



ac 



(42) 



replaces the definition the same quantity in eqs.([l4|). The Hamiltonian op- 
erators are given by 



f-\D 



J'r 



\Da 



Db 



\ 



\cD + c x 



\aD \{Db + bD) 
\ bD \Dc-c x {b 2 - ac)D + D(b 2 - ac) ) 

D 3 \ 

—D 2 a D 

D 2 bD + DbD 2 
+DaDaD j 

and the corresponding densities of Hamiltonian functions 

1 





( 





J'l = 





D 3 






D 3 


-DaD 2 



H' = c. 



a{D- l b) 2 - (D- l b) (D^c) 



(43) 



(44) 



(45) 



yield the system ( |I~4D with this redefinition of e. These Hamiltonian operators 
are also compatible. The symplectic 2-form obtained from the inverse of (|4|) 
is given by 

J = dp A dr - q xx dp A dp x + p xx dp A dq x + ^dq A dq x (46) 
which in terms of / in eq.(|6|) is simply 



u' = 2 f tt df tx A dfu + — df x A df tx 



(47) 



where we have again discarded a total derivative. 
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12 Covariant formulation 



We have noted at the beginning that eqs. @ and @ are obtained by a simple 
flip of the independent variables. However, any comparison of the Hamilto- 
nian operators (||^) and (f44|) with fllTf ) and (|18|) yields nothing more than 
a complete mismatch even though they arise from what is in fact the same 
WDVV equation. The fact that the Hamiltonian operators for these equa- 
tions look very different forces us to look for a unifying framework which 
exists in the covariant formulation of symplectic structure [IIU] . Indeed these 



different looking Hamiltonian structures are simply different components of 
the Witten-Zuckerman closed, conserved current 2-form. In order to show 
that this is indeed the case we start with the Lagrangian for, say eq.(||), 

= Ijftt ftxx ~ ftx fttt ~ ~^fxx ftx (48) 

expressed in terms of the original variable /. This is the same as the La- 
grangian ( |26|) up to a total derivative. The Witten-Zuckerman current 2-form 
a/' which is closed and conserved 

0/^ + 0^ = 0, d(J* = 0, n = x,t (49) 



can be obtained directly from the Lagrangian (48). We find 

uj x = 2 f tt df tx A df tt + \df x f\ dftx, 
w* = 2 f tx df tt A df tx — \ df x A df xx 



(50) 



which satisfies the properties of the symplectic 2-form listed in eqs.(§^). If 
we go back to eqs.(|58]), P7| ) for the symplectic 2-forms appropriate to eqs.(^), 
(||) expressed in terms of / alone, then we find simply 

uj x = u', u 1 = u, (51) 

as we should expect, since these WDVV equations are related by a flip of t and 
x. Unfortunately the lack of a variational formulation of WDVV equations 
d§|), or (^) that yields their first Hamiltonian structure through Dirac's theory 
of constraints makes it impossible to present a similar covariant description 
of their first symplectic structure. 
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13 Second pair of WDVV equations 

In order to discuss the Hamiltonian structure of the WDVV equation (|TT1) 
we need to redefine 

b 2 + c 2 + be — ac — a 

(52) 



e 



a ■ h 



in the system of evolution equations (|i~4"|). The bi-Hamiltonian structure of 
these equations of hydrodynamic type is obtained through the same process 
we presented earlier and the result can be summarized in the following. 

Theorem 2: 

The equations of hydrodynamic type (|14]) with e given by (H) admit 
bi-Hamiltonian structure ([IB]) with Hamiltonian operators 



Jo 



I -aD - Da 

\{a-b)D-Db \{ 

V bD -De 



\D(a -b)-bD Db - cD \ 

>-c)D + \D{b-c) \(c-e)D + Dc 
\D(c-e)+cD 



lD + De + eD ) 



Ji 



-D 3 
D 3 



—D 3 D 



D 3 
—D 3 

a + b + 1 jj2 
a + b 



—D 3 
D 2 a + b + l D 
a + b 



\ 



D 



~a+T) 



D b+C +^D 



V 



-D 



a + b 

D b + C M + ^ D-^D 
a + b a 

a + b + 1 jj a + 



a + b 



tu 



D 



(53) 



(54) 



and the Hamiltonian functions which are integrals of the densities Tio = c 
and 

U x = (c-b-a) D^b D^c+ib+l) D^a D- 1 b+cD- 1 aD- 1 c+-{D- 1 a) 2 (55) 

2 



yield the equations of motion. The Hamiltonian operators Jo and J\ are 
compatible so that the theorem of Magri is applicable in this well. 

Comparison of these results with the bi-Hamiltonian structure of eq.([10|) 
presented in again results in a mismatch which, however, is not as severe 
as before since eqs. (|T0|) and flTT|) resemble each other rather closely. In order 
to show that these different Hamiltonian structures simply correspond to 
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different components of the Witten-Zuckerman 2-form first we note that the 
variational principle with the Lagrangian 

C= ~ \ (f x l + ft. /xx) + fxxx ftt fxx + ftxx ftx ftt (5g) 
fxxx ftt ftx fttt ftx f xx 2 fttt ftx 

yields eq . (|TTD , or more precisely we get a linear combination of both t and 
x derivatives of eq.flTlf). This Lagrangian can be rewritten in terms of the 
auxiliary variables (|24]) 

£ = (-Qx Qxx + Qx r x -r p xx + p x q xx + \p x ^j p t - (§g£ + p x q x ) r t 

+ (q x r x - q x p xx + Px Pxx +p x r x ) q t - q x r r x - \p% (57) 

~^~Qx Qxx V Qx Pxx T Px Qx Qxx Px V fx Qx Px 

which yields a triplet of evolution equations, the integrability conditions of 
which result in eq.(|TT[). 

The symplectic 2-form obtained from the inverse of the Hamiltonian op- 



erator (54) is given by 



& = l(r x + q xx ) dp x + (r x - p xx ) dq x - (p xx + q xx ) dr } A (dp + dq) 
— \dp A dp x 

and this is rather similar to the result 

= {(r x + q xx ) dp x + (r x - p xx ) dq x - (1 + p xx + q xx ) dr ] A (dp + dq) 
+ \dq A dq x 

(59) 

for eq.fllOD that was obtained earlier in 0. The Witten-Zuckerman symplec- 
tic 2-form that follows from the Lagrangian (|56D is 

u x = 2 (j txx - f ttx - |) df tx A df x - (f ttx + 1) df xx A df x 
+ (2 ftxx - fxxx) dfu A df x + (2 f tt + df xx ) df xx A df tt 
+2 ftt dftx A dfu - fttx dftx A df t + f txx df tt A df t - 2 f tt df tx A df xx 
= u' 

to 1 = \df xx A df x + [(fttx + ftxx) df xx + (fttx — fxxx) dftx 
- (fxxx + ftxx) dfu] A (dfx + dft) 
= uj 

(60) 
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where to, Co' are given by eqs.(|59"D and (^3) respectively. This is exactly the 
same final result flSTD we found earlier for eqs.(§) and @. 

14 Conclusion 

The bi-Hamiltonian structure of WDVV equations of associativity is ex- 
pressed by different pairs of Hamiltonian operators depending on which in- 
dependent variable in the free energy is chosen to play the role of time. 
However, we have shown that the symplectic 2-forms which can be obtained 
from the second Hamiltonian operators are simply different components of 
the covariant Witten-Zuckerman conserved current symplectic 2-form. We 
have been able to obtain this result only for the second Hamiltonian structure 
because only in this case does there exist a degenerate Lagrangian subject to 
second class constraints which through the Dirac bracket yields the Hamilto- 
nian operator. The lack of a variational principle for the WDVV equations 
that yields their different first Hamiltonian structures makes it impossible to 
discuss them together in the framework of the covariant theory of symplectic 
structure. 
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